Formation kinetics of magnetic chains, rings, X's, and Y's 
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The kinetics of aggregation in a monolayer of magnetic particles are studied using stochastic 
dynamics computer simulations. At low densities (< 8% coverage) the equilibrium structure is 
made up of chains and rings; the primary mechanisms by which these motifs form are described. 
At higher densities (> 15% coverage), we observe large transient concentrations of Y-shaped and 
X-shaped defects that ultimately give way to an extended, labyrinthine network. Our results suggest 
that a defect mechanism - such as that proposed by Tlusty and Safran, Science 290, 1328 (2000) - 
could drive a metastable phase separation in two dimensions. 
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The structure, phase behavior, and dynamics of 
strongly interacting, dipolar fluids present considerable 
challenges to soft-matter physics. The most common re- 
alization of a dipolar fluid is a ferromagnetic colloidal sus- 
pension, or ferrofluid. In the ideal case, this consists of 
spherical, homogeneously magnetized monodisperse par- 
ticles with diameters ~ 10 nm, sterically stabilized and 
immersed in a nonpolar solvent. The resulting colloidal 
interactions are caricatured by the widely studied dipo- 
lar hard sphere fluid. Despite almost four decades of 
intensive experimental, theoretical, and simulation study 
at least one outstanding question remains to be an- 
swered definitively: are point dipolar interactions alone 
sufficient to drive vapor-liquid phase separation? 

On the one hand, the Boltzmann-weighted, angle av- 
erage of the dipolc-dipole potential gives (to leading or- 
der) an isotropic, attractive pair potential that varies like 
— 1/r^, where r is the interparticle separation; this is 
expected to produce conventional condensation behav- 
ior p|. On the other hand, simulations show that con- 
ventional condensation is preempted by strong aggrega- 
tion, driven at low temperatures by the energetically fa- 
vorable 'nose-to-tail' conformation Q. If phase separa- 
tion occurs in 3D, then it is of a rather unusual nature; 
simulations suggest that the low-density phase mainly 
consists of chain-like aggregates, while the high-density 
phase resembles a network of interconnected segments |5| . 
One possible scenario involves a defect-mediated phase 
transition in which the chains' defects are the singly- 
connected particles at the chain ends, while the network's 
defects consist of particles with three or four near neigh- 
bors in a Y-shaped or X-shaped conformation, respec- 
tively. Such fundamental issues are not only of relevance 
to magnetic fiuids; the physical properties of many ma- 
terials are governed by connectivity and topology j^]. 

Recently, 2D dipolar fluids (with 3D magnetostatics) 
have received attention due to the possibility of directly 
imaging aggregation in thin films using cryogenic trans- 
mission electron microscopy 0, 1 ^ llfll . The equilibrium 
structure 0, 0| and dynamics |l3l | of 2D dipolar flu- 
ids have been studied in detail using computer simula- 



tions. At low density and low temperature, the domi- 
nant structural motifs are isolated chains and rings. The 
high-density structure consists of a labyrinthine network 
of long chains, with a small concentration of X-shaped 
and Y-shaped defects. There has never been any sugges- 
tion of a vapor-liquid phase transition in the 2D system, 
but a transition between isolated and system-spanning 
clusters at low-density has recently been characterized 
[l^ . The relative ease with which complex self-assembled 
structural motifs can be visualized and analyzed in exper- 
iments and simulations means that 2D fluids are attrac- 
tive and important objects of study. Given the signif- 
icance of such systems, it is surprising that more isn't 
known about the self-assembly process itself, starting 
from a 'random' arrangement of particles. In one of 
the few studies in this area. Wen et al. directly imaged 
the aggregation of nickel-plated glass microspheres 1^. 
They observed that rings can be formed by two short 
chains joining at both ends simultaneously. Branching 
at X and Y defects was not considered. 

In this Letter, we present a detailed simulation study 
of the aggregation process in monolayers of strongly dipo- 
lar particles. Starting from equilibrated configurations of 
non-polar particles, we elucidate the mechanisms of clus- 
ter formation that occur when the dipoles are 'switched 
on'. At low densities, rings are predominantly formed 
by single, isolated chains folding up (rather than by the 
association of two short chains). Interestingly, at high 
densities, the system shows a high transient concentra- 
tion of defect particles. This is significant because it sug- 
gests that highly branched structures could be stabilized 
kinetically. It might therefore be possible to realize a 
metastable, defect-mediated phase separation in the lab- 
oratory. We also briefly consider the types of mechanisms 
by which the structure evolves at long times. 

We model the system as a monolayer of monodisperse 
dipolar soft spheres. The interparticle potential is given 
by 
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where e is an energy parameter, cr is the sphere diameter. 
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fii is the dipole vector on particle i, r is the interparti- 
cle separation vector, and r = \r\. Reduced units are 
defined as follows: temperature T* — k^T/e, where fee 
is Boltzmann's constant; dipole moment ^* = \J /i^/ ta^] 
particle number density p* = pa^ , where p ~ N/L"^ and 
L is the length of the square simulation cell; time t* = i/r 
where r — -^/mo^/e is the basic unit of time. Stochastic 
dynamics simulations were performed according to the 
integrated Langevin equations [l^ 
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where Fi [T.i) is the net force (torque) acting on dipole i 
at time t, fii = fJ-t/p is a unit dipole orientation vec- 
tor, St is the integration time step, and Dq {Dq) is 
the translational (rotational) diffusion constant at infi- 
nite dilution. The components of the 2D vector SWl 
and the 3D vector SW[ were generated independently 
from Gaussian distributions subject to the conditions 
(SWl) = 0, {SWl ■ SW^) = ADl6tS,j, {SW[) = 0, and 
{SW[ ■ SWj) — GOgStSij. In this scheme, the short-time 
inertial dynamics are suppressed: with high dipole mo- 
ments these occur on timescales of order 1 (in reduced 
units) ll."^ . while of primary interest here are conforma- 
tional processes occurring on timescales orders of magni- 
tude longer. 

We present results for TV = 1024 particles at a range 
of densities, with a large dipole moment p* — 2.75, 
and temperature T* — 1. Characteristic diffusion con- 
stants were estimated from the (stick) Stokes-Einstein 
laws yielding = kBT/S-Krja ~ 4 x 10~^^ m^ s~^ and 
Dq ~ k-QT/Tirjcr ~ 1 x 10^ s^^ for spherical particles 
with (7 = 10 nm, in a solvent of viscosity = 10~^ Pa s 
at temperature T = 300 K; the dimensionless quanti- 
ties D^r/o-^ 0.004 and D^t = 0.01 were obtained 
using the mass for 10 nm spheres with mass density 
~ 8000 kg m~'^ (typical for iron or cobalt) and energy 
parameter e = k^T. The integration time step was 
5t* — 0.01. Self-assembly was initiated from configu- 
rations generated with p* — 0. For each density studied, 
five independent runs with different initial configurations 
were conducted, and the results for each density were 
averaged. The configurational temperatures associated 
with the positions and orientations of the particles were 
measured independently In ^ill cases the instan- 

taneous configurational temperatures fluctuated about 
T* = 1, with rms deviations of about 0.1, throughout 
the self-assembly process. 

At equilibrium, the low-density structure mainly con- 
sists of small rings and chains |l2j. In Figs, ^a)-(d) 
we show how a ring is formed by a chain closing in on 
itself, at a density p* = 0.05. Wen et al. suggest that 




FIG. 1: (Color online) Instantaneous configurations at var- 
ious densities and times: (a) p* = 0.05, t* = 3500; (b) 
p* = 0.05, t* = 4330; (c) p* = 0.05, f = 4800; (d) p* = 0.05, 
t* = 4950; (e) p" = 0.5, t* = 10000; (f) p* = 0.5, t* = 50; (g) 
p* = 0.2, t* = 1080; (h) p* = 0.2, t* = 1200; (i) p* = 0.2, 
t* = 1300; (j) p* = 0.2, t* = 1400. The black lines between 
particles denote 'bonds' identified using an energy criterion. 
In (a)-(d) and (g)-(j) the orientations of the dipole moments 
are indicated with red and blue cones. 



rings are formed by two short chain-like segments making 
connections at either end simultaneously 15]. In con- 
trast, movies of the aggregation process in our simula- 
tions show that the majority of rings are formed from 
isolated chains. At high density the equilibrium struc- 
ture resembles a labyrinth of long, winding chains |l3l| : 
an 'equilibrium' configuration at p* — 0.5 is shown in 
Fig. n^e) . Locally, the equilibrium structures at low and 
high density are not so different; most particles are in 
chain-like environments, flanked by two near neighbors. 

The process of self-assembly was monitored by identi- 
fying particles belonging to the same cluster on the basis 
of an energy criterion. We chose an energy cut-off of 
Uc = — 0.6/i*^, which captures a range of likely confor- 
mations for neighboring dipoles, either in chains or in 
more exotic defect environments We identified 'ter- 
minal particles' (one near neighbor), 'internal particles' 
(two near neighbors), 'defect particles' (three or more 
near neighbors), chains, rings, and 'defect clusters' (con- 



FIG. 2: Numbers of clusters containing (from top to bottom) 
n — 1, 2, 4, 6, 8, or 10 particles, as functions of time: (a) 
p* = 0.1; (b) p* = 0.2; (c) p* = 0.3; (d) p* = 0.5. The solid 
lines are simulation results, and the dashed lines are fits to 
the von Smoluchowski model. 

taining at least one defect particle). 

In Fig.|21we plot the numbers of n-mers (n = 1, 2, 4, 6, 
8, and 10) as functions of time at densities in the range 
0.1 < p* < 0.5 ll^. To analyze the simulation results, 
we use a simple von Smoluchowski model |l8| 

^ijP^P] - Pn ^^niPi (4) 
i+'j—n i—1 

where p„ is the number density of n-mers, and the fcy's 
are 'rate constants'. Under the simplifying assumption 
that all rate constants are equal (kij = k) , the number of 
n-mers is 7V„ = 7V(M)"~V(1 + kptY+'^, and the total 
number of clusters is Y^^=i = + kpt), where p 

is the total number density of particles. We determined 
a reduced rate constant k* = kr/a'^ for each density by 
fitting to the total number of clusters, with the results 
k* = 0.049(3), 0.104(5), 0.156(4), 0.247(6), and 0.52(1) 
at p* = 0.05, 0.1, 0.2, 0.3, and 0.5, respectively. The 
resulting curves for Nn{t) are shown in Fig.|3 The gen- 
eral level of agreement is quite good at short times, but 
at longer times the von Smoluchowski model underesti- 
mates the numbers of clusters. This is likely to be due to 
complex and slow relaxational processes occurring as the 
structure ripens, and the gross simplification of setting 
all rate constants equal. 

In Fig. 121 we show the numbers of chains, rings, and 
defect clusters as functions of time at densities in the 
range 0.1 < p* < 0.5. At all densities, the numbers of 
chains are greater than the numbers of rings and defects 
during the early stages of the aggregation process, but 
the long-time behaviors are very different. Rings are fa- 
vored over defect clusters at low density (p* < 0.1), in 
accord with the known equilibrium structures. At high 
density {p* > 0.3) defect clusters are favored over rings; 



FIG. 3: Numbers of chains (excluding monomers) (dia- 
monds) , rings (circles) and defect clusters (triangles) , as func- 
tions of time, at densities (a) p* = 0.1, (b) p* — 0.2, (c) 
p* = 0.3, and (d) p* = 0.5. 



the numbers of defect clusters show strong maxima at in- 
termediate times {t* ~ 10-100). At intermediate density 
(p* — 0.2) there are roughly equal numbers of rings and 
defect clusters. 

To further elucidate the aggregation mechanisms, in 
Fig. ^ we plot Xn{t), the fraction of particles with n 
neighbors {n =^ 0-4) at time t, at densities in the range 
0.1 < p* < 0.5. (No particles with five or more near 
neighbors were observed.) At all densities, xq falls mono- 
tonically, xi shows a maximum, and a;2 — > 1 as the ma- 
jority of particles ultimately end up with two neighbors 
as parts of chains or rings. Dramatic differences between 
low-density and high-density aggregation kinetics are ev- 
idenced by X3 and X4. At low density {p* < 0.1) less 
than ^ 0.2% of particles have three or more neighbors 
at any given time. At intermediate density (p* = 0.2) 
X3 shows a maximum at t* ~ 100, while X4 is essentially 
negligible. At higher densities (p* — 0.3, 0.5) and at in- 
termediate times (t* 100), up to 10% of all particles 
are defects, with approximately ten times more particles 
having three neighbors than four neighbors. In Fig. ^f ) 
we show a snapshot of a configuration at p* — 0.5 and 
t* = 50, i.e., close to the location of the peak in 2:3. Note 
the higher connectivities of particles within clusters as 
compared to those when the system is closer to equilib- 
rium [t* = 10000, Fig. me)]. Ultimately, the number of 
defect particles falls again by at least an order of magni- 
tude as equilibrium is approached. 

Of course, aggregates continue to disassemble and 
form, even at equilibrium. During the aggregation pro- 
cess at high density, there is a very slow net decrease in 
the number of defect particles within the system. It is 
impossible to identify specific local events where the num- 
ber of defects is reduced irreversibly, but in Figs.nig)-(j) 
we show a sequence of snapshots that illustrate the kinds 
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FIG. 4: Fractions of particles with n neighbors, as functions 
of time, at densities (a) p* = 0.1, (b) p* = 0.2, (c) p* — 0.3, 
and (d) p* = 0.5: n = (solid line); n — 1 (circles); n = 2 
(diamonds); n = 3 (triangles); n — 4 (squares). 

of processes which, in the long run, lead to a changing 
cluster distribution. The time interval t* = 1080-1400 
corresponds to where xs is decreasing relatively rapidly. 
In the upper regions of these figures, we see how two 
chains segments meet, bond temporarily, and then sepa- 
rate. In the lower halves of Figs.^^g) and^h), we see 
how two chains can collide at right angles, and split in 
to two new chains. In Fig.^i), we see a cluster with two 
defect particles; in Fig.^j) we see the result of that clus- 
ter having ruptured at the positions of those two defects 
almost simultaneously. 

In summary, our results show that the aggregation 
mechanisms at low density and high density are very dif- 
ferent. At low density, aggregation proceeds through the 
formation of chains, some of which will go on to form 
rings, as shown in Fig. n^a)-(d). At high density, it ap- 
pears that the transient aggregates include a significant 
proportion of defect clusters, such as those implicated in 
3D dipolar phase separation 0]. At equilibrium, the par- 
ticle connectivities at low density and high density are 
somewhat similar in that the vast majority of particles 
have exactly two neighbors. 

Our observations suggest a direct experimental test for 
the theory of defect-mediated phase transitions put for- 
ward by Tlusty and Safran . If the ferrocoUoid particles 
were modified to introduce a strong, short-range attrac- 
tion (e.g. chemically, or with added polymer to induce 
depletion forces) then the transient network structure at 
high density may be kinetically stabilized long enough 
for phase separation to occur, even in thin films where no 
equilibrium transition is anticipated. This can only arise 
if the high-density branched structure corresponds to a 



local free-energy minimum. The experimental approach 
suggested here could provide a way of guiding the dense 
phase in to that minimum, but there may be alternative 
strategies. In any case, we hope that aggregation kinetics 
and metastable phases in strongly dipolar fluids can be 
investigated by direct experimental observation. 
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